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Abstract 



The Navier-Stokes-Fourier model for a 3D thermoconducting viscous fluid, 
^" 
^H ' where the evolution equation for the temperature T contains a term propor- 

_ ^ ' tional to the rate of energy dissipation, is investigated analitically at the light 

'c/2 ' 
^^' of the rotational invariance property. Two cases are considered: the Cou- 

Oh' 

ette flow and a flow with a radial velocity between two rotating impermeable 

►^ and porous coaxial cylinders, respectively. In both cases, we show the exis- 

m 

ff^ ' tence of a maximum value of T, Tmax, when the difference of temperature 

o' 

r^ ' AT = T2 — T1 on the surfaces of the cylinders is assigned. The role of Tmax is 



discussed in the context of different physical situations. 
FAGS numbers: 47.32.Gc 

In Ref. PI, a thermoconducting incompressible viscous fluid system, named Navier- 



K^ ! Stokes-Fourier (NSF) model, is presented. In 3D, this model is governed by the equations 

Ut + u-Vu + V(-)-z/V2u = 0, (1) 

P 

V ■ u = 0, (2) 

Tt + u-VT- knV^T = -^J^i^^uj + d.u.f, (3) 

where u is the velocity fleld, T = T{x, y, z, t) is a passive scalar identifled by the fluid tem- 
perature [@], p = p{x,y, z,t) the fluid pressure, p the fluid density, r] the dynamic viscosity. 



V = rj/p the kinematic viscosity, Cp the heat capacity, k, the heat conductivity, kn = K,/{pCp) 
the thermal diffusivity, and di = d/dxi {xi = x, X2 = y, x^ = z). 

The r.h.s. in (^) is related to the rate of energy dissipation e = 2uS'^, where S"^ = SijSij 
and Sij = {l/2){diUj + djUi) are the strain matrix elements. Among the many questions 
arising in the study of the model (|I|)-(|^), such as for instance the investigation of statistical 
solutions and the onset of turbulence [^ , the presence of a dissipation term in the equation for 
the temperature deserves a special attention. At the best of our knowledge, in three spatial 
dimensions Eqs. (HI)-® have never been studied by exploiting an analytic procedure. In 
this Letter we discuss some new thermal effects described by exact solutions of the vortex 
type of the NSF model. In particular, we have analyzed a steady Taylor-Couette flow for 
an incompressible fluid entrapped between two rotating coaxial impermeable cylinders, and 
a flow with an additional radial component of the velocity, when the walls of the cylinders 
are porous (for technical details on the experimental devices, see @,|^). As is well-known, 
the hydrodynamic problems connected with the system of two coaxial rotating cylinders has 
been intensively studied from an experimental and theoretical point of view since the work of 
Couette (1890) |^ . The main motivation was the discovery by Taylor of unstable flow regimes 
consisting in toroidal vortices (1923), which appear for high values of the Taylor number [|^. 
This number can be deflned in several forms. Following ||5| we will use Ta = riQid/u, where 
Vi is the radius of the inner cylinder, d the width of the gap between the two cylinders, Qi 
the rotational velocity of the inner cylinder. In our context, we will consider low Taylor 
numbers {Ta ~ 100 -r- 200), in order to preserve the geometrical symmetries of the model. 

In the two cases of impermeable and porous cylinders, we show the existence of a value 
^max for the temperature fleld associated with the fluid, when the difference of temperature 
AT = T2 — T1 on the surfaces of the cylinders is assigned for asymptotic values of time. This 
maximum appears as a consequence of a mechanism of energy dissipation in heat, due to 
viscous effects, and depends on the width of the gap between the cylinders, on their relative 
angular velocity, and on the flxed value of AT. Remarkably, we found that the position of 
^max as a function of the radial coordinate r can vary continuously in the gap ri < r < r2 



between the two cylinders by varying the parameter AT. 

In the NSF model, p is taken to be constant. Although this assumption could appear 
as a too severe restriction on the possible modelling of realistic situations, nevertheless the 
NSF model deserves to be analyzed as a guide-framework for the investigation of less ideal 
systems, such as the Rayleigh and the Lorentz models, which describe incompressible fluid 
motions where thermal phenomena are more significant ||l[|. 

We found that the NSF model allows Lie-point symmetries generated by infinitesimal 
operators giving finite group transformations and exact solutions via the corresponding re- 
duced equations. (The technical machinery yielding these reductions is outlined in @J^). 
Many of these symmetries have a nice geometric interpretation: they express the invariance 
of the NSF system under rotations, translations, Galilean boosts and scale transformations. 
Therefore, as one expects, their expressions do not explicitly contain the viscosity v. Never- 
theless, such symmetries hold only in the viscous case, in the sense that the corresponding 
symmetries for the inviscid situation cannot be reproduced by those for z/ 7^ by setting 
z/ = 0. 

Here we shall limit ourselves to deal with the rotational symmetry defined by the gener- 
ator 

Vr = ydx - xdy + U2du, - Mi<9„2 . (4) 

A more systematic treatment of the symmetry properties of the NSF model together with 
the investigation of the role of the boundary conditions will be reported in a separated paper. 
The operator (H) gives rise to a reduced NSF (RNSF) system whose last equation is 



9t + Ui9r + U39, - -^ { 2Ut, + 2-i- + Ut + Ui 

„2 < -^2z^^Zr^ ^^■iz^'^^lzVZr 



'^'^'' + % + Ul + Ul + 2Ul + 2UM 






kn Ur + J + 9zz] = 0, (5) 



where z, t, r = ^Jx^ + y^ are independent symmetry variables and \J\, U2, t/3, n, 0, related 
to the original variables by 



Ui = Ml COS (f + U2 sin (f, U2 = —Ml sin ip + U2 cos ip, 



(6) 



are dependent variables expressed in terms of z, t, r. Here cosip = x/r, simp = y/r, Ui, U2 
can be interpreted as the radial and the azimuthal components of the (reduced) velocity U, 
and f/3 is the component along the 2;-axis. 

We find that the RNSF system admits two solutions which are, respctively, i) the circular 
Couette flow (in which only the azimuthal component f/2 of U is different from zero), and 
a) a flow where all the three components of U are different from zero. In this case t/s turns 
out to be time dependent in such a way that f/3 — > exponentially when the adimensional 
variable ^ot goes to infinity (see ([T^)). In the original velocity components wi, U2 and M3 
the circular character of the Couette flow appears, while the flow of case ii) shows a spiral 
behavior. 

Case i) 

A simple but notable solution to the RNSF system is obtained by choosing [/i = f/3 = 
and 



U2 = - + (3r, 
r 



(7) 



where a,/3 are constants. Then 



P 



2r2 2 



+ -/3V^ + 2a/3plnr, 



(8) 



apart from an arbitrary function of time, while the temperature T is provided by 



T = exp(-Aot) X 



ci-^o I \ —, -r I + C2Y0 



'Ao — Ai 



k 



H 



Ai , . Ai 
C3 cos \ - — z + C4 sm \ - — z 
\ kH \ ku 



V kn 



X 



-d) 



(9) 



where 



is the contribution to the temperature coming from the energy dissipation present in (^ 
and Ao > Ai > 0, A2, Tq and Cj are constants. Here Jo and Iq are the Bessel functions of 
the first and the second kind, respectively (see JlOl, P- 358). 



The original velocity components corresponding to the reduced velocity U = (0, -+/3r, 0) 



are 



Ui 



-(^ + /3)y, n2=(^ + /3)x, «3 = 0, (11) 



from which we have r = tq = const. Equations (|Tl]) are solved by x = tq cos Kt, y = 
Tq sin Kt, z = const where K = a/r^ + (3. 

The reduced solution (|^ is pertinent to a steady flow where the fluid occupies the gap 
ri < r < r2 between two coaxial cylinders of radii ri and r2 rotating with (constant) angular 
velocities Qi and ^^2 [|ll|]. Under the hypothesis of no-slip condition on the cylinders, a 
and P take the form a = [r^r|(fi2 — fii)]/(r^ — r|), /? = (r^fii — r|fi2)/('^i — 1^2)^ while the 
vorticity 0;^^ associated with the (reduced) velocity (|^ coincides with the vorticity uj for the 
velocity (|lT]). 

In Fig. |l| we show the behavior of the temperature field, as being described by the 
solution d^), as a function of t and r (brighter regions correspond to higher values of T). By 
exploiting the scale invariance of the NSF system we can express the physical constants in 



an adimensional form |T2|. Hereafter we will assume z/ = 10 , p = I, r] = 10 ,C„ = 1, 



K = 6x 10^'^. After a transient phase, only the contribution ([I0|), due to the dissipative term 
in (^) survives. In this situation we have a maximum of T for any t in the gap between the 
two cylinders. This maximum, as shown in Fig. ^ for the asymptotic case, can continuously 
vary by varying AT (or equivalently A2). 



FIGURES 




FIG. 1. Temperature field for tlie motion of a fluid between two coaxial rigid and impermeable 
cylinders with radii ri = 1 and r2 = 2 for z = 0.5 (q = 1). 
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FIG. 2. Asymptotic behavior of the thermal field for l<r<2, 0<z<5, when AT = 
{Ta = 150). 



In Fig. ^ we plot T — Ti vs. r for different values of 6 = ^, in the case AT = 0. We 
point out that when the cylinders rotate in the opposite sense (6 < 0), T^ax increases. This 
is related to a growth of the energy dissipation induced by friction effects. 



I-Tl 



b=-0.1 
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FIG. 3. Behavior of the thermal field as a function of the radial distance between the two rigid 



and impermeable cylinders for some values of the parameter & = q7 with ri = 1 and r2 = 2 
{Ta = 150). 



Case ii) 

Another interesting vortex solution of the NSF equations (0)-(0) can be found by putting 



into the RNSF system, where A, B and 7 are constants. 
The pressure p = II is given by 

p(4A5r^+2 - 752 - 7^2) pA2r2(^+i) 



(12) 



P 



+ 



27r2 ' 2(7 + 1) ' 

apart from an arbitrary function of time, while the RNSF system yields for f/3: 

f/3 = exp(— o-Q^)^^ X 
^0 J- ixl^r] + A,Y. (Mr " 

, 2 vv z/ / 2 vv z/ 



(13) 



(14) 



where ao > 0, Aq, Ai are constants. 

An interesting physical interpretation of the vortex solution presented in Case ii) is the 
following. Let us consider a device constituted by two rotating porous coaxial cylinders of 
radii ri < r2 containing in the gap a fluid of kinematic viscosity u. 

Let 7 = ^^^^^^ be the radial Reynolds number, and mq = Ui{ri) the radial velocity through 
the wall of the inner cylinder M. We assume that the flow is inward for 7 < and outward 



for 7 > 0. For Aq = Ai = 0, the velocity filed reproduces the solution investigated in [Q. 
Furthermore, for 7 = the quantities (|T2|) correspond to the Couette flow (0). In the case of 
porous cylinders, A and B are related to the geometry and the dynamics of the device, and 
explicitly become: A = -[ilia'^il - b/a'^)]/[r^{l - a^+^)] and B = rl^iil - ba')/{l - a'+^), 
where a = ri/r2 and b = ^72/^1- 

In general, i.e. for f/3 7^ 0, the evolution equation (|^) for the temperature is very 
complicated and the finding of exact solutions is a difficult task. However, for Aq = Ai = 0, 



we obtain 



2ku 




, r + C2Y^nL. 




Cssm 



//iO - yUl 



k 



Z \ + C4 COS 



/yUo - yUl 



H 



A; 



// 






(15) 



where 



II 



2,,2\ 



T. 



/i2r'=ff 



+ 



Ar]AB 



r"< + 



2r]{B^ + 7^ 
pCp{2kH + Ji^)r^ ' lpCp{kH - // 

2pCp(7 + l)[2A;^(7 + l)-7'^] 






(16) 



Here Hq> jii> 0, /i2 and Tq^^ are constants, the Prandtl number Pr = u jku is supposed to 
be 7^ 1 and 7 > (i.e. the radial flow is considered outward) to avoid singularities. 

To determine the time behavior of the coordinates (a;(t), y{t), z{t)) of the fluid particle, we 
integrate the velocity fleld u = {ui, U2, Ms) to get r^ = 2Ct + rg, where x and y obey the pair 
of equations of the time-dependent oscillator type which afford the solutions x = y/rcosip, 
y = ^/Tsm'^p where r = 2Ct + r^ and ip = ^[^^2+^ j^- B\n^]. 

The term (|16D is again due to the presence of the dissipation rate in (^ , and represents 
the asymptotic limit of the solution (p!5D, for t -^ 00. In Fig. ^ we plot the temperature 
fleld (|I5]) for z = 0.5. This plot is analogous to that of Fig. (|1|). In other words, also in this 
case a maximum of T emerges after a transient phase. 

The position at which Tmax is located in the gap ri < r < r2 (Fig. |^) depends on the 
value of 7, on the difference of temperature AT and quadratically on the width of the gap. 



In particular, Tmax tends to migrate towards the well of the outer cylinder when the values 
of 7 increase and vice-versa. A similar situation occurs when AT changes. 

In Figs. ^ and |^ we plot the thermal field for b = ^ > and b < 0, when AT = and 
AT = 4. 




FIG. 4. As in Fig. || for the case of porous cylinders (7 = 1) 
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FIG. 5. Asymptotic behavior of the temperature field between two rotating porous cylinders 
as a function of r (1 < r < 2) and z (0 < z < 5) for AT = 0, 7 = 5 and 6 = {|^ = 0.2 {Ta = 100). 
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FIG. 6. Temperature behavior as a function of the radial distance between the two rotating 
porous cyhnders for several values of 6 = ^ and 7 = 1 if the cylinders temperatures are kept equal 
{Ta = 100). 
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FIG. 7. As in Fig. |, where now AT = 4 and 7 = 1 {Ta = 200). 

To conclude, from our analysis of the temperature behavior of vortices of the NSF model 
(|l|)-(0), notable thermal effects arise as a consequence of the presence of the energy dis- 
sipation term in the evolution equation (^ for the temperature field. We hope that our 
theoretical results can stimulate some experimental work addressed to a possible evidence 
of such effects. 
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